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We analyze the phenomenology of a prolonged early epoch of matter domination by an 
unstable but very long-lived massive particle. This new matter domination era can help 

o ■ 

to relax some of the requirements on the primordial inflation. Its main effect is the huge 
entropy production produced by the decays of such particle that can dilute any possible 

> : 

unwanted relic, as the gravitino in supersymmetric models, and thus relax the constraints 



on the inflationary reheating temperature. A natural candidate for such heavy, long-lived 
particle already present in the Standard Model of the electroweak interactions would be a 
heavy right-handed neutrino. In this case, we show that its decays can also generate the 
' observed baryon asymmetry with right-handed neutrino masses well above the bound from 

^ , gravitino overproduction. 

I. INTRODUCTION 
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■ of the Universe, that were not in causal contact in the radiation last scattering era, have approxi- 
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mately the same temperature today [2]?. In inflationary models, an epoch of exponential expansion 
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Inflation was introduced in the 80s [1] as a solution to several problems of the big bang cosmol- 
ogy. Perhaps the main problem was the large-scale smoothness problem. Why different patches 



inflates a small patch of the Universe in causal contact to contain all the observable Universe to- 



day. Simultaneously if the temperature after inflation is low enough, inflation helps also to dilute 
unwanted relics from higher scales and reduces the flatness problem. 

It is usually assumed that some kind of inflation starts already at the Planck scale to avoid 
the Universe collapse in a few Planck times if £1 > 1 or (for any Q) to prevent the invasion of 
the surrounding inhomogeneity to our homogeneous patch before inflation. On the other hand, 
the scales observable today in the cosmic microwave background left the horizon at an energy 
U 1 / 4 < 6 x 10 16 GeV , or 60 e-foldings before the end of inflation. So that, inflation must 
end below this scale. After the end of inflation comes an era of reheating when the inflaton field 
oscillates around its minimum and decays to ordinary particles. The final reheating temperature 
where we recover ordinary big bang cosmology can take any value from U 1 / 4 above to scales as low 
as 1 MeV. However, the required dilution of unwanted relics, as GUT monopoles or gravitinos in 
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supersymmetric models, forces the reheating temperature to be well below the GUT scale or even 
below T RH < 10 8 GeV in SUSY models. 

In this letter, we propose a simple and economic mechanism that helps solving some of these 
problems and reduces the requirements on the primordial inflationary mechanism without further 
additions to the particle spectrum of the Standard Model with right-handed neutrinos . After 
an initial inflationary epoch (still necessary to reproduce the observed correlation on temperature 
fluctuations at large scales) we assume our Universe is radiation dominated for a short period and 
then enters a matter domination era due to the existence of a heavy long-lived unstable particle 
that decays to radiation well-before nucleosynthesis, when we connect with usual cosmology. In 
the Standard Model, as we will show, this role could be played by a heavy right-handed neutrino. 
In the literature, it is well-known that late time entropy release can help to ameliorate some of the 
problems of standard cosmology. However, so far most of these works have only considered moduli 
fields in supersymmetric theories (see for instance (J, 0, f], ^Hsl) and their real presence in nature 
could be considered more speculative than the existence of right-handed neutrinos. 

As we show below, this matter domination mechanism, naturally embedded in the SM, is able 
to help primordial inflation in several aspects. A long period of matter domination can reduce 
mildly the number of e-folds before the end of inflation at which observable perturbations were 
generated, relaxing this way flatness conditions on the inflationary potential. Moreover, the large 
entropy production in the decay of this particle completely dilutes any unwanted relics, eliminating 
the constraint on the inflation reheating temperature. In this sense our matter domination epoch 
has the same advantages as thermal inflation [9], without resorting to yet another scalar field and 
/or scalar potential. 

II. REWRITING THE HISTORY OF THE UNIVERSE 

Let's assume for a moment that at an early time a massive particle dominated the energy 
balance of the Universe by many orders of magnitude. How does the observed Universe feel this 
new epoch?. Are there an y o bservable consequences of this? 



As it is well-known 



id, 



a massive particle X becomes non-relativistic when the temperature 
of the thermal bath falls bellow its mass, Mx, and its energy density freezes out when it drops 
out of equilibrium. Then, Y-relic abundance relative to photons (radiation) becomes constant 



1 From now on, we call the Standard Model with right-handed neutrinos simply "Standard Model" 
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and therefore, if X is completely stable, the energy density of the X particles will eventually 
become larger than the radiation one, dominating the energy balance of the Universe. If, X is 
not completely stable but rather sufficiently long-lived to dominate the energy density, later on X 
will decay into relativistic particles that thermalize. The radiation content of the Universe will be 
increased and it will enter again in a radiation dominated era. This will be the basic evolution of 
the Universe in our model. 

We start from a radiation dominated Universe at a scale larger than the mass of our X-particle, 
Mx, where this particle is in thermal equilibrium. We assume that this particle decouples from 
the plasma at a temperature of the order of its mass. This particle is unstable, although very 
long-lived, i.e. it has very weak interactions with radiation degrees of freedom. Then, its energy 
density, px, starts diluting as matter, much slower than radiation. If its lifetime, Tx, is long 
enough, it will necessarily dominate the energy density of the Universe. Although our X particle 
decays all the time through an exponential law 12(], it is only when the age of the Universe is of 



the order of 1/Tx that the decay will sizeably reduce its abundance. Once it reaches this point, 
it will fastly decay into radiation and our Universe will go back to a radiation dominated epoch 
where it will connect with the usual cosmology. This "matching" with the standard scenario must 
happen well before nucleosynthesis. 

The evolution equations for the matter and radiation energy-densities are well-known: 

px = -3H(l + w x )px ~ T x px (1) 
pf d = -AHpf d (2) 

ft™ = -AHp^ + TxPx (3) 
2 d 8tt f old | „new , „ A (A \ 
H ~ a ~ M^S Pr Pr { ) 

where wx is the equation of state parameter of particle X, which drops from 1/3 to zero as X 
becomes nonrelativistic, p° ld is the energy density in radiation not related to X decays, p]} ew is 
the energy density in radiation produced by X decays and here we assume a flat Universe after 
inflation consistent with observations 

From these equations we can see that, when X is nonrelativistic, the number of X's per comoving 
volume (Nx = R 3 px/Mx) follows a simple exponential decay law and the (formal) solution to 
these equations is given by: 

-3 



PX = p x 



aP 



e -r x (t-t o) (5) 
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FIG. 1: Evolution of the different components of the energy density of the Universe from T ~ 10 GeV 
to T ~ 1 MeV. Short-dashed (green) line corresponds to pj?^ versus time, long-dashed (blue) line to p new 
and the solid (red) line to p x . We start from p x = p° ld /200 at T ~ 10 15 GeV, with a T x = 1CT 20 GeV. 
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where the superscript zero denotes the value of that quantity at the initial epoch. 



In general, it is not possible to integrate analytically these equations, although useful approx- 
imations exist [rJ . However it is always possible to solve them numerically as we have done 



to generate Figures [T] and [2j These e vo 
M. Turner and collaborators in Refs. 



ution equations, ([T]) to were thoroughly analyzed by 
12l . [r3 | (for a more recent work see for example ) and we 
agree completely with their analysis of the matter and radiation densities, entropy and tempera- 
ture. However, we are specially interested in the particular limit of small Tx and large Mx- In 
fact, as we will see below, we will focus on the limit of small Tx keeping Tx ^> 1/issiv s ° that no 
trace of X is present at nucleosynthesis time, in agreement with observations. 

In this limit, our massive particle X dominates the energy density of the Universe by many orders 
of magnitude during a sizeable fraction of the thermal history of the Universe (see Figure [T]) . As 



shown in Refs. 



12 



13], during the X decays the temperature of the Universe does not fall as i -1 / 2 
(a^ 1 ), but rather as t -1 / 4 (a~ 3 / 8 ) due to the entropy release of the decays and the temperature 
reaches an almost flat plateau from the point where the energy density in new radiation born 



5 




FIG. 2: Temperature of the Universe (in GeVs) as a function of time. The four different epochs in the 
evolution of the Universe, radiation domination, matter domination, decay and radiation domination again 
can be seen in the different slopes of the curve. The time dependence of the temperature is explicitly 
indicated for each epoch. 

through X decays and the energy density in old radiation become comparable up to t ~ T^ 1 . 
After this time, X rapidly decays and the temperature falls again as t -1 / 2 (a -1 ) (see Figure [2]). 
Once it has completely decayed away, our Universe is left with a temperature 



Wdecay = 1.0-10' [— ) GeV (9) 

where g* counts the effective number of relativistic degrees of freedom. Notice that the temperature 
after A-decay depends only on the decay width T. The ratio of entropy per comoving volume before 
and after X decay is given by 

Spre-decay \ T J V 10 GeV / 

where r = gx/% if X is a boson and r = 3 gx/8 if it is a fermion, with gx the total number of spin 
degrees of freedom. 

As mentioned before, an early period of matter domination, triggered by a long lived massive 
particle that goes out of equilibrium and comes to dominate the energy density of the Universe 
before decaying, reduces the number of e-foldings before the end of inflation at which our present 
Hubble scale equaled the Hubble scale during inflation, i.e. the time of horizon crossing. The 
reduction is given by 

AiV = lln( Ppost - mattcr - dN ) (11) 
12 \ Ppre-matter-d / 
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FIG. 3: Comoving Hubble radius, log (1/aH), versus log a. This plot shows the different eras entering the 
e-foldings calculation. Inflation is an epoch where log(l/ai?) is decreasing. Exponential inflation gives 
a line with a slope of -1. In all other cases the inflation line is shallower. During matter domination 
(1/aH) oc a 1 / 2 , while during radiation domination (1/aH) oc a. The current dark energy domination 
signals a new inflationary epoch. The horizontal (black) solid line indicates the present horizon scale. The 
number of e-foldings before the end of inflation at which observable perturbations were born is the horizontal 
distance between the time when (1/aH) first crosses that value and the end of inflation. The solid (red) line 
represents a Universe with a period of matter domination before BBN. The dashed (blue) line represents 
the standard cosmological history of the Universe with only one (recent) epoch of matter domination. 

where /9p re -matter-d an< ^ /°post-matter-d are t ne energy densities at the beginning and end of the mat- 
ter dominated era, respectively. The reduction can also be expressed in terms of the X parameters 
as 

^-imf^'V^) (12) 

This reduction is illustrated in Figure O where it can be clearly seen that to determine the number 
of e-foldings after horizon crossing of a given cosmological scale, as the present Hubble scale, the 
complete thermal history of the Universe must be used. From nucleosynthesis onwards this history 
is well in place. However earlier epochs are still very uncertain. The standard cosmological model 
assumes that inflation gives way to a long period of radiation domination (we neglect here the 
period of reheating that immediately follows inflation and assume sudden transitions between the 
different regimes). The radiation dominated epoch lasts until a redshift of a couple of thousands 
before entering an era of matter domination, which at redshift below one gives way to the current 
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acceleration. 

Changing the sequence of events after inflation can therefore have a strong impact on the the 
number of e-foldings calculation. If our Universe goes through a long period of a regime where 
log(l/aH) scales as a n , i.e. H oc a~~( n+l \ it is straightforward to see that with n > 1 the total 
number of e-foldings will be increased while for n < 1 this number will be reduced. A period of 
matter domination belongs to the latter class, as in a matter dominated epoch (1/aH) oc a 1 / 2 
opening the door to a significant reduction on the number of e-foldings. 

If we put all this information together we can see what are the required features of our particle 
X, its mass Mx and its lifetime T, if it is to dominate over the energy density of the Universe for 
a long period. 

Although we do want a prolonged period of matter domination, we want it to come to an 
end at the latest shortly before nucleosynthesis, as the Universe must have attained thermalized 
radiation domination by that time. Using Eq. ([9]), this condition sets a lower bound on T, 

1 /2 

T > 2.0 • 10- 24 (|0 GcV . (13) 

We can also get un upper bound on T, by requesting the reheating temperature Tp 0S ^ eca y to be 
at most 10 8 GeV, so that no unwanted relics will be produced after X decay in supersymmetric 
models. Such a condition reads 

T < 0.7 ( —) 1 GeV . (14) 

V2ooy V ' 

However, for lifetimes this long, we can see from Eq. (|10p that only large X-masses are capable of 
effectively diluting the unwanted relics. In the case of the gravitino this even more difficult as the 
gravitino abundance is also proportional to the temperature, T > Mx- 

Of course detailed bounds can be set only after specifying the basic physics behind X, its 
production mechanism, its decays, or in short, its interactions. Nevertheless, we can say that 
requiring at least five orders of magnitude of dilution by entropy production for Mx < 10 10 GeV 
would need 

/ a \ 1/2 / a \ 1/2 

2.0 • 10~ 6 — GeV > T > 2.0 • 10~ 24 — GeV . (15) 
V200,/ ~ ~ V200/ v ' 

In this case the five orders of magnitude of increase in the entropy should be enough to get rid 
of unwanted relic which could have been produced at earlier times for ? 1 p 0St _ c j eca y < 10 8 GeVs. 
Larger lifetimes are also possible if we do not need such a large entropy production. 
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With regards to the reduction in the number of e-foldings, only a very prolonged period of 
matter domination, i.e. large Mx and V in the lower part of the allowed range, is required to give 
a significant reduction. In most cases, however, this number is expected to be below 10. 

III. MATTER DOMINATION IN THE STANDARD MODEL 



Now, we must check whether this early matter domination epoch could exist in the context of 
the Standard Model of the strong and electroweak interactions or any of its extensions. Clearly to 
obtain such an early matter domination era we need a massive particle, X, in thermal equilibrium 
at temperatures above its mass with a very long lifetime. The simplest candidate for our X-particle 
would be a right-handed neutrino. Right-handed neutrinos are one of the minimal additions to 
the Standard Model to reproduce the observed neutrino masses through the seesaw mechanism 
[15I . These right-handed neutrinos, Ri, have super-heavy masses, that can be as high as the Grand 
Unification scale, and are singlets under the SM gauge group. The only renormalizable couplings 
of Ri with the SM particles are, possibly small, Yukawa couplings. Therefore, if these couplings 
are sufficiently small, it seems possible that the right-handed neutrinos play the role of X-particle 
with a long lifetime. 

More precisely, the right-handed neutrino masses and Yukawa couplings have to reproduce the 
measured neutrino masses and mixings though the seesaw mechanism, m UL = v \ Y v ■ (Mr) -1 • Y„ , 
with V2 the vacuum expectation value of the up-type Higgs. From here it is straightforward to 
obtain the required right-handed Majorana matrix from the seesaw formula itself: 



Mr = v\ Yj ■(rn VL y X -Y v 



(16) 



From the light neutrino mass matrix, m UL , we know the mixings and the two mass differences. The 
mixing matrix U is close to the so-called tribimaximal mixing, 

/ * 73 \ 
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V y/E V3 V2 ' 



Then, m 
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U* • Diag (mi, m-2, 7713) • U\ and the inverse of this matrix is 
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Diag ( ■ U T . Therefore, 
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(17) 
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(18) 
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Experimentally we have that 7713 = m atm ~ 0.05 eV, 777.2 = m so \ ~ 0.008 eV and mi <C m-2 in 
the normal hierarchy situation and 7772 = m a t m — 0.05 eV, mi = m atm — m so \/2 ~ 0.046 eV and 
7773 <C m,2 in the inverse hierarchy case [lfj]. 

As seen in Eq. (|16p the masses of the right-handed neutrinos reproducing the observed light- 
neutrino masses and mixings are determined by the Yukawa matrix, Y v . Choosing the basis of 
diagonal Y^Y V and diagonal charged lepton Yukawa matrix, we have Y v = Vl ■ Diag (2/1, 2/2, 2/3). 
Obviously the physics depends strongly on the form of Y u , both on the eigenvalues, yi, and the 
Vl matrix. Let us first analyze the role of Vl- We have two limiting situations: a) Vl has large 
mixings and is the source of the observed PMNS matrix in neutrino mixings, Vl — U* and b) the 
mixings in Vl are small, similarly to the situation observed in the CKM mixing matrix, Vl — 1. 

Case a) is very simple, the seesaw mechanism plays no role in the generation of the neutrino 
mixings. The observed large neutrino mixings are already present in the Yukawa couplings before 
the seesaw mechanism. This corresponds to the situation where, the light neutrino Majorana mass 
matrix and the neutrino Yukawa couplings, or equivalently, the right-handed neutrino Majorana 
matrix and the Yukawa combination Y*Y, can be simultaneously diagonalized. So, we have Mr = 
v\ Diag (7/1 /mi, 7/2/7772, 2/3/7713). The decay widths of the right-handed neutrinos will be given by 
Tj = g^rMj ^Y^Y^) and we must compare it with the Hubble rate, H(T = Mi) in order to know 
if/when our massive neutrino will go out of equilibrium. Equivalently, we can compare the effective 
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mass 77ij = (y^Y^J v 2 /Mi (i.e. r = ^ ^f) with the critical mass m» = lx 10~ 3 eV 
right-handed neutrino would dominate the energy density if rhi < m^/g^ where g* is the number of 
radiation degrees of freedom at T = M$, The presence of g* is due to the fact that in a time H^ 1 
we can see from Eq. ([5]) that the ratio of matter and radiation densities grows as a. But matter 
has to overcome g* radiation degrees of freedom and hence to dominate the energy density it needs 
a lifetime (/* times longer 2 . Altogether, in case a), it is clear 7713 = 7773, 777,2 = m2 and ttii = mi. 

Therefore, in the normal hierarchy case, taking mi < 10 -10 eV, Ri would dominate the energy 
density of the Universe with a mass of, 

This case is therefore a perfect example of how a right-handed neutrino can dominate the energy 
density of the Universe after inflation. In terms of rhi we can write Eq. (|10p as 

Vst-dccay ^ Q 4 r ( M ^ ( 1 x IP' 6 oV \ ^ 
Spre-decay ^ 200 ^ V mi J 



However, as we see below, this large effective mass, rhi — 10 5 eV, does not generate sufficient entropy. 
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Therefore, if we want two orders of magnitude of entropy production, we would need mi = 10 10 
eV (c_ di n g to P ^4. 6 X UT» GeV) a»d r p08 , decay a 2 X U^Q.V. Here, a, the gtavitino 
abundance is approximately linear with the reheating temperature [19(, two orders of magnitude 
of dilution by the entropy production would correspondingly relax the bound on the reheating 
temperature by two orders of magnitude. Naturally, this can be easily improved by choosing 
smaller mi and y\ in Eqs. (|19p and (I20j) . 

However this simple model has several phenomenological problems. First, given that fh\ <C m*, 
this right-handed neutrino would not be produced thermally through its Yukawa interactions. In 
fact this will be a common problem of any massive particle dominating the energy density of 
the universe as necessarily its decay /production rate will be much slower than the Hubble rate. 
Therefore we will always need another active interaction to produce our right-handed neutrino 
in the thermal plasma after inflation. This role could be played, for instance, by a gauged B— 
L interaction. Many Grand Unified models based on 50(10) or groups containing it have an 
intermediate scale of the order of 10 13 GeVs with a intermediate gauge group containing U(1)b-l, 
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221 ] . In these grand unified models 



as SU(2) L x SU{2) R x SU{3) C x U{1) B -l for example 
the B-L coupling unifies with the other gauge couplings at Mqjjt and therefore it is always strong 
enough to keep the right-handed neutrinos in thermal equilibrium in the unbroken phase. However, 
the B-L gauge interaction can never mediate the neutrino decay as it couples only diagonally in 
flavour. The neutrino decay would require a Yukawa interaction to lighter states that, as shown 
above, in this case is very small. A second problem, more specific of this particular case is that 
the decay of this right-handed neutrino erases completely any previously existing baryon or lepton 
asymmetry and therefore, we need some mechanism to generate the observed baryon asymmetry. 
The baryon asymmetry generated by this completely out-of-equilibrium decay of the right-handed 



neutrino is given by 



8 1 x „ 

^ = 23^ 16^ E 

37^ 



Im 




Mi 


( 







(21) 



but, as Y^Y is completely diagonal in the basis of diagonal right-handed neutrino masses, no new 
lepton asymmetry is generated by Ri decays. 

However, this situation is very unstable and a slight departure from the perfect alignment of 
the Majorana and Yukawa matrices changes the situation. If we call R the rotation diagonalizing 
the neutrino Yukawas in the basis of diagonal left-handed neutrino Majorana masses, we have 
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Y u = U* ■ R ■ Diag (2/1,2/2,2/3) and, 

1 „ „ 1 „ „ 1 



(Mr) „ = «| 2/iJ/j flii-Ry + i? 2i i? 2i + «3i^3j • (22) 

'mi 771 2 7T13 



If the Yukawa couplings are sufficiently hierarchical, 2/3 ^ 2/2 3> 2/1 , the heaviest eigenvalue will be 
given by the element (Mr) 33 , 

(M r ) 33 ~ t; 2 y 3 2 f— (sin0i 3 ) 2 + — (sin0 23 ) 2 + — (cos i3 cos 23 ) 2 ) , (23) 
\mi m 2 777,3 / 



24]. From here, we can see that 



where we used the standard PDG parametrization for the matrix R 
the contribution from m\ « m 2 ,m3, will dominate (Mr) 33 , and hence the heaviest right-handed 
neutrino eigenvalue, if (sin0i3) 2 > 7711/7713. For m\ = 1 x 10~ 10 eV, a sin #13 > 4.4 x 10~ 5 will be 
enough and such a small departure from perfect alignment will completely change the situation. 
To understand this, it is enough to analyze a simpler situation with 0i 2 = 2 3 = and #13 7^ 0. Let 
us take 7/3 — 1, 2/1 — 10~ 6 (similar to the up-quark hierarchy), 7713 ~ 0.05 eV and ra\ ~ 10~ 10 eV. 
In this case, the two right-handed neutrino eigenvalues (the other one is unchanged) are given by, 

R<f 22/ cos 2 % sin 2 6>i3 \ 

M R3 - V 2V3 + , 

\ rn 3 mi J 

M Rl ~ v\y\ T —— 

mi cos #13 + )7i3 sin t/13 

(24) 

If sin 0i3 3> 7711/7713 then both M3 and 7713 are fixed by mi, while M\ and mi are fixed by 7713. So, 
we have, for mi in the interesting range: 

7773 = mi/ sin 2 0i3 (25) 

And this is the only neutrino that can dominate the energy density of the Universe. Clearly, this 
situation is not interesting phenomenologically as it is not possible to produce it after inflation and 
(probably) it does not produce a large amount of entropy. 

The most interesting situation corresponds to sin0i3 < mi/m 3 . In this case from Eq. (|24p 
we have Mr z ~ v\ y^/ms (l + sin 2 #13 7713/7711) and Mr 1 ~ v\ y\j (mi (1 + sin 2 0i3 777,3/777,1)) . 
Now, the rotation that diagonalizes the right-handed mass matrix is given by sin<^> ~ 
sin 0i3 (2/17773) / (2/377*1) . Therefore in the basis of diagonal right-handed neutrino masses we have, 

/ y\ (\ + sin 2 0i3 (0) 2 ) -2/12/3 Si sin 13 \ 



yty 



(26) 



V -2/12/3 ^ sin 0i3 2/| / 
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and for sin #13 < 7711/7773 we have rh\ ~ m\ and 7713 ~ 7713. This means the lightest right-handed 
neutrino, with a mass approximately given by Eq. (I19p . can still dominate the energy density. In 
such a situation, we can see from Eq. (I2ip . that the generated baryon asymmetry is given by 

2" 



Im 

1 

rjB 



771-3 

2/12/3 — sin Wi3 



7711 ~ J J yi m3 ~ 1 „,2 oi „2 fl m 3 



U»7T / y 2\ mi 167T m ^ 



^ sin^ia -§. (27) 



Taking 7/1 ~ 1CP 7 , 7773 ~ 0.05 eV, mi ~ 10~ 12 eV and sin #13 = 0.1 x 7771/7773, we would obtain 
t]b — 10 _7 sin99 with tp the CP violating phase of (Y'Y)^. Therefore, in this case, it would be 
possible to generate the observed baryon asymmetry and simultaneously dilute the relic density 
and in particular the gravitino density by three orders of magnitude. This means that the bound 
on the inflationary reheating temperature would be relaxed by three orders of magnitude. Clearly, 
using smaller values for mi and y\ the situation can be improved arbitrarily. 

The second example of Yukawa mixing matrix was case b) where Vl — 1 so that we can neglect 
this rotation on m UL as a small rotation will not modify the contribution of the different neutrino 
eigenvalues to the matrix elements. Then, we have: 



to,, 1 — m, 1 m 1 —m, 1 



M R = v 



2 



2 m o 1 +2m. 1 ,,. 

( 2/1 3 2/12/2 3 2/12/3 3 ^ 

mZ 1 —m7 1 o rnT 1 +2mZ 1 +3mZ 1 mr 1 +2mr 1 -3m7 1 

2/12/2 3 2/2— 1 1 2/22/3 — 1 *~ 

rn~ 1 —m7 1 mT +2m^ 1 —3m7 9 mr 1 +2mT 1 +3m^r 1 , 

V 2/12/3 3 2/22/3 — 1 — V3— 1 — / 



(28) 



3 6 «3 6 

We must diagonalize this matrix to obtain the right-handed neutrino eigenvalues and the Yukawa 

matrix in the basis of diagonal Mr. In analogy with the charged lepton and quark Yukawas we 

can expect 2/3 S> 2/2 ^ 2/1- Then we obtain, in the normal hierarchy case, 

Mr 3 ~ v 2 yl ^j- , M R2 ~ t; 2 yf 27773 1 and M «i - « 2 l/l 3 ™2 ( 2 9) 

Then, we have, 7773 ~ 6 mi, 7772 — 7773 and 7771 ~ | 7772. This means that once again the right- 
handed neutrino that could dominate the energy density of the Universe is the heaviest one and its 
mass would be given by Eq. ([25]) with sin #13 ~ l/v6, i-e. a mass close or even above the Planck 
scale, unless 2/3 is much smaller than 1. 

Thus we see that with hierarchical Yukawa eigenvalues, similar to the up-quark eigenvalues, it 
is possible to have right-handed neutrino dominance of the energy density with consistent phe- 
nomenology, although only in rather fine-tuned situations where Vl is very close to the PNMS 
mixing matrix but not exactly equal. If we move to an extension of the minimal Standard Model 
with right-handed neutrinos (or Minimal Supersymmetric Standard Model) the same results can 
be obtained without such a tight fine-tuning. 
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From Eq. (|16p . we can see that besides the Yukawa mixing Vl we can still use different Yukawa 
eigenvalues. In a Grand Unified Theory (GUT) with an underlying Pati-Salam symmetry, we 
would expect the neutrino Yukawa couplings to be related to the up quark Yukawas, and in fact 
we would expect one of the neutrino eigenvalues of the order of the top Yukawa coupling 251 ] . 



However the two light Yukawa eigenvalues are less restricted. The masses of the right-handed 
neutrinos will depend on the Yukawa eigenvalues and we can make them as small as we wish. 
However, if y% is large, we will normally be in the same situation as before and the only neutrino 
with a sufficiently large lifetime will still be vr 3 , which will be far too heavy. An interesting limit 
is when one of the Yukawa eigenvalues is exactly zero, y\ = 0, and therefore, one of the light left- 
handed neutrino masses is zero. Again the simplest situation is when the right-handed Majorana 
matrix and the Yukawas, YJY V are simultaneously diagonalizable. In this case, it is clear that only 
two right-handed neutrinos will play a role in the seesaw mechanism and the third one will be 
completely decoupled from the seesaw. In fact this third neutrino does not couple to the doublets 
through Higgs Yukawa couplings. Given that right-handed neutrinos are singlets under the SM 
group, apparently these neutrinos do not decay at all. However, if we have a GUT symmetry, as 
for instance SO(W) or a group containing SU(2)r, at a high scale, the right-handed neutrino will 
decay with a lifetime, 

r^^f-, W(^) 5 (^C^) 4 ,30, 

So, indeed we can see that this right-handed neutrino can dominate the energy density of the 
Universe if it is produced through another interaction, as a gauged B-L. In this case the production 
of the baryon asymmetry is not possible through gauge interactions. However, it is possible that 
this right-handed neutrino has non-vanishing complex Yukawa couplings when it unifies with the 
quarks at the GUT scale 3 . In this case the neutrino decay through GUT Higgses could violate CP 
and generate the observed baryon asymmetry. 

Finally we would like to point out another possibility where we introduce more SM singlets 
mixed with the three "standard" right-handed neutrinos. An example of this situation is provided 
by the so-called "double seesaw" [27[ mechanism. In this case, the right-handed Majorana masses 
are generated through a second seesaw with these additional singlets. The singlets would decay 
only through its mixings with the three right-handed neutrinos and, as we are introducing another 
free parameter, we can easily make any of these new singlets to dominate the energy density of the 



3 For instance, we could think of a Georgi-Jarlskog vev distinguishing up-quark and neutrino Yukawas and being 
zero for the neutrinos [261 ] . 
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universe. Similarly, these singlets would easily generate the observed baryon asymmetry through 
the usual neutrino Yukawa couplings. The problem of how to generate a thermal abundance of 
these singlets could be solved again if they are charged under a gauged B-L symmetry. 



IV. CONCLUSIONS 



In this work, we have shown that an early epoch of matter domination by a long-lived massive 
particle can help to solve some of the problems of primordial inflation. We have seen that the large 
entropy production generated by the decays of such particle can dilute unwanted relics from higher 
temperatures, relaxing the constraints on the inflationary reheating temperature. In supersym- 
metric theories this mechanism can help to solve the gravitino problem. Moreover, a long period 
of matter domination reduces the number of e-foldings before the end of inflation at which the 
observable cosmological perturbations were generated. In the Standard Model a natural candidate 
for such heavy, long-lived particle is a heavy right-handed neutrino. For low enough mass of the 
lightest left-handed neutrino and neutrino Yukawa mixings sufficiently close to the PNMS mixing 
matrix, the right-handed neutrino dominates the energy density of the universe for a long time and 
generates a large amount of entropy in its decay. In this case, we show that its decays can also 
generate the observed baryon asymmetry for right-handed neutrino masses well above the bound 
from gravitino overproduction. 
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